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ABSTRACT 

We study tangent bifurcation of band edge plane waves in nonlinear Hamiltonian lattices. 
The lattice is translationally invariant. We argue for the breaking of per mutational symme- 
try by the new bifurcated periodic orbits. The case of two coupled oscillators is considered 
as an example for the perturbation analysis, where the symmetry breaking can be traced 
using Poincare maps. Next we consider a lattice and derive the dependence of the bifur- 
cation energy on the parameters of the Hamiltonian function in the limit of large system 
sizes. A necessary condition for the occurence of the bifurcation is the repelling of the band 
edge plane wave from the linear spectrum with increasing energy. We conclude that the 
bifurcated orbits will consequently exponentially localize in the configurational space. 
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I. INTRODUCTION 



An increasing number of publications deal with the phenomenon of discrete breathers 
in nonlinear translationally invariant Hamiltonian lattices (for a history on the subject see 
|l||, 0, 0, [^], 0, 0, 0, p and references therein). Discrete breathers are time-periodic 
solutions of the equations of motion. These solutions are exponentially localized in the con- 
figurational space (see and in detail ||TU|)- One can view discrete breathers as analoga of 
the breather solutions of the sine-Gordon (sG) partial differential equation |^. Breathers of 
the sG equation are structurally unstable (nongeneric) against perturbations of the Hamilto- 
nian density of the sG field ||l2l. In contrast discrete breathers appear to be structurally 
stable (generic) solutions, with the seeming neccessary ingridients being the nonlinearity and 
the discreteness (i.e. the existence of a finite upper bound of the linear spectrum) |13|. One 
could think that the sG PDE is the 'bottleneck' which gives contact to the discrete breathers 
by simply discretising the spatial differential operator. This thought is most probably wrong, 
because in order to maintain structural stability of discrete breathers one has to avoid reso- 



nances of multiples of the breather's frequency with the linear spectrum [Ty]. Thus discrete 
breathers appear to be rather unrelated to their counterpart in the sG case. 

Proofs of existence of discrete breathers have been given so far for i) arrays of weakly 
coupled anharmonic oscillators by use of analytical continuation of periodic orbits |T^ and 
for ii) chains with homogeneous interaction potentials by use of discrete map analysis of the 
Fourier coefficients (with respect to time) flSl . Structural stability of discrete breathers has 
been shown for one-dimensional lattices by use of discrete map analysis of Fourier coefficients 

HI- 

It turns out to be rather complicated to give a proof of existence for discrete breathers for 
an arbitrarily choosen Hamiltonian function. Still for applications or numerical studies it is 
important to know whether a given system can posess discrete breathers. From the Fourier 



coefficient analysis in [jTO[ we can predict, that if the frequency of the breather comes close 
to the linear spectrum, then the exponent of the spatial decay becomes very small. If the 
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energy of the breather solution becomes small in the same limit, then the breather solution 
comes close to the solution of a band edge plane wave (here band edge plane wave stands 
for a plane wave with a wave vector of the band edge of the band of allowed frequencies of 
the linearized equations of motion). Consequently one could expect that discrete breathers 
appear due to tangent bifurcations of band edge plane waves if one increases the energy Ezb 
of the plane wave starting from Ezb = 0. Indeed the birth of discrete breather solutions 
as bifurcating periodic orbits due to the tangent bifurcation of band edge waves has been 
observed numerically for one- dimensional Fermi-Pasta-Ulam (FPU) lattices ||15| with the 



number of lattice sites ranging between = 40 and = 100. Note that tangent 
bifurcation implies that a pair of Floquet multipliers of the linear stability analysis of a 
periodic orbit collide at +1 on the unit circle. 

The bifurcation problem of band edge plane waves has been a problem studied inde- 
pendently of the discrete breather phenomenon Budinsky and Bountis [|16| have 



approached the problem using Floquet theory. They have then analyzed the eigenvalue prob- 
lem for large N. However there were discrepancies between the results in |TB| and the study 
by Sandusky and Page |T5|. Also in both cases only Fermi-Pasta-Ulam (FPU) chains were 
considered. Further no analytical results about the properties of the new bifurcated orbits 
were given. Consequently the goal of the present paper is i) to derive exact conditions for 
the tangent bifurcation of band edge plane waves for large A^, ii) to consider both FPU and 
nonlinear Klein-Gordon (KG) systems, and iii) to show that the bifurcated orbits will break 
the permutational symmetry of the system and exponentially localize in the configurational 
space. 

In section II we analyze the bifurcational problem for a system of two coupled oscillators 
(A^ = 2). With the help of Poincare maps we will show that the derived formula apply. 
In section III the general lattice with A^ degrees of freedom is introduced, and the periodic 
orbits of the band edge plane waves are obtained. The results of the bifurcational analysis 
are presented in section IIIA,B. In section IIIC we show that the bifurcated orbits break the 
permutational symmetry. We discuss our results in section IV. 



II. THE CASE OF TWO COUPLED OSCILLATORS 



Let us consider the following Hamilton function 

H = ]^Xl + ]^Xl + V{X^) + V{X2) + ^{X^-X2) , (2.1) 
which describes the dynamics of two coupled oscillators using the equations of motion 

The potential functions from ( |2.1| ) can be expanded in the form 

OO 



OO -| 

^J,=2 



(2.2) 
(2.3) 



Note that we have to demand f 2 > and 02 > in order to ensure the potential energy in 
(p.l|) being a positive definite quadratic form in the limit of small energies. 



A. The normal mode periodic orbits 



For small energies we can neglect the anharmonic terms /i > 2 in ( |2.2| ),( P^ . In that 
case it is appropriate to transform the original variables into normal coordinates 

n = ^(Xi + X2) , = ^(Xi - X2) . (2.4) 



Using the relations given above we obtain the equations of motion for the normal modes 
including all anharmonic terms of the potential functions (|2.2|) ,( ^73D : 



Yo 



00 (i^<(M-i)) 

E 

^=2 u=l,3,... 



00 (!^<(M-1)) 

E^M E 

^=2 u=0,2,... 



( 



/i 



V 



2 ' 



V 



11=2 



(2.5) 



(2.6) 
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Note that the sums over v in these equations run over even or odd values of v respectively 
because cancellation of terms due to (p.4[). 

We can now solve ( p .51 , 12 .61 ) immediately for the two periodic orbits, which correspond to 
the standard harmonic normal mode solutions in the limit of small energies: 

oo 

I: F2 = o , yi = -^t;,rr' , (2.7) 

[1=2 

oo 

II: yi=0 , % = -Y^v;^r^ , t;2..+i = . (2.8) 

^=2 



Here the renormalized constants for case II in ( p.8|) are given by 



v^ = v^^ 2^- • (2.9) 



Note that we had to require f2m+i = in case II ( p.8|) , since in the case of nonvanishing odd 
terms in the potential Vi^z) the right hand side of (|2.5| ) has terms containing the variable I2 
only, and the solution would be characterized by nonzero Y\. 

The solutions of the differential equations ( p.7| ),( pl8|) are given in terms of elliptic func- 



tions. Here we will be interested in the limit of small energies (amplitudes) and apply 

standard perturbation theory ( [|T8| Chapter 7). For case I (equation ( |2.7| )) we obtain 

+00 



A;=— 00 



^ 2t>2 3^2 

Ao = -^Ei , A\ = ^E^ , A2 = -p^i?i , (2.12) 



f2 2t>2 6f; 



2 

For case II (equation ( |2.8| )) we obtain in a similar way 



2 

Ex = Iy,^ + V{Yi) . (2.13) 



hoo 



fc= — oo 



1^2 = E ^fce*'"* ' = A.k , (2.14) 
2^^ ~ 3^1 



u^^ = v.+ ('^-'S]E. , (2.15) 



A, = -^E2 , Aj = -^E, , A2 = ^E2 , (2.16) 

V2 2V2 6^2 

^2 = -Yi + ViY^) , Viz) = E . (2.17) 

^ AJ=2 



All terms not present in ( |2.10|) -( ]2.17|) are of order or higher and can be neglected, as 
we will see in the following. 



B. Tangent bifurcation for periodic orbit I 

Let us consider a small perturbation 61^2 of the periodic orbit I 

^1,2 ~^ + ^1,2 • 

If we linearize the resulting equations of motion for the perturbation we obtain 

00 

5i = -Y.ili-l)v,Yr'Sr , (2.18) 
^=2 

00 

^2 = - E(/^ - ^)^^^Yr^^2 - 2</)2<52 . (2.19) 

ti=2 



Note, that the variable Yi in ( ^.18[ ),( P?TUD is the time-periodic solution of (^.7]). As we see. 



the perturbations do not couple with each other. Equation ( |2.18| ) describes the continuation 
of the periodic orbit ( |2.7| ) along itself (shift of time origin) or along the one-parameter family 
(change of energy Ei or frequency uj). The associated pair of Floquet multipliers is obviously 
located at +1 on the unit circle [Ol. 



Equation ( p.l9|) then describes the nontrivial perturbations (the associated phase space 
is two-dimensional, and can be visualized with the help of Poincare maps). Let us rewrite 

^2 = -[(t^2 + 202) + /(t)]52 , (2.20) 

00 

fit) = Y.(^^-l)v,Yr' . (2.21) 



This type of equation is called Hill's equation [|T9|, since the function f(t) is periodic in time 
(with the period of Yi): 

f{t) = f{t+-) . 
We are interested in tangent bifurcations which appear when 



y/V2 + 202 



(2.22) 



The function f(t) in ( |2.21| ) can be expanded in a Fourier series: 

+00 

fit) = Y: F,e''^' , Fik) = F.k . 



(2.23) 



If the function f{t) becomes infinitesimally small, then according to the Floquet theory 
parametric resonance can appear for the cases 



"2" ' ^ • 



(2.24) 



With ( p.22| ) the relevant pairs (A;, k') become (1, 2) , (2, 1). In the limit of small values of 
El we can then retain in ( p.20|) the significant terms only: 



5o = - 



+ 2Ficos(cjt) + 2F2Cos(2cut) 



=V2 + 



(2.25) 
(2.26) 



Using relations (p.ll| - p.l3|) we obtain for small Ei 



V4 



-2^ + 3— El 



V2 V2 
V3 



Fi = V2 
\3vi 2v2 



(2.27) 
(2.28) 

(2.29) 



A tangent bifurcation will appear if the solution for the perturbation 62{t) becomes periodic 
with the period of the function f{t) (since at the bifurcation the associated pair of Floquet 
multipliers has to merge at +1 on the unit circle): 

For given potential functions ( p.2| ),( p^ this can happen only at certain energies E^. Ap- 



plying the method of strained parameters |18| to our problem we obtain two solutions 



i) : 02 = , El arbitrary , 



9f2f4 



(2.30) 
(2.31) 



Since we have to require positive values for Ei and 02 tangent bifurcation can take place 
only if 

f4 10 fo , , 

- < TT^ 2.32 

f 2 9 f 2 

is fulfilled. Remarkably this condition ( p.32|) is equivalent to requiring the frequency of 



the periodic orbit I being a decreasing function of Ei (cf. ( |2.11|) ). Note that we have 



to exclude the nongeneric case of f 3 = f 4 = (however we have no restrictions on 03^4,... 
since the anharmonic interaction terms simply do not contribute to the linearized equations 
(^),(g39D). The two solutions (p:30|)-(p3lD define two lines in the space {^1,02}- The 



area between the two lines corresponds to the parameter cases when the periodic orbit is 
hyperbolic due to tangent bifurcation. Note that the requirement of small energies implies 
also small values of 02- For large values of Ei and 02 higher order corrections will appear. 

C. Tangent bifurcation for periodic orbit II 

In analogy to the previous chapter we derive the linearized equations for the perturbation 
of the periodic orbit II: 

00 

51 = -Y.if^-l)v^Yr'5^ , (2.33) 

fi=2 
00 

52 = -Y.if^-l)v,Yr'52 . (2.34) 

fj,=2 

Note that we had to require V2m+i = 0. Again the perturbations do not couple with each 
other. Equation ( |2.34]) describes continuation of the periodic orbit II. 



Equation ( |2.33| ) gives the nontrivial perturbations. It is again a Hill's equation. In the 
limit of small energies E2 we obtain in analogy to ( p.20| )-( p.29| ) 

S^ = - + 2F2Cos(2cut)] Si , (2.35) 
n^ = V2 + Fo , (2.36) 
Fo = 3^E2 , (2.37) 

V2 

F2 = . (2.38) 

2f2 
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Requiring periodicity Si(t) = Si(t + 211/00) and applying the method of strained parameters 
we find the two solutions for a tangent bifurcation 

■ ^2 = ^^.2 L ^ > (2-39) 
8003 ~ 36r;204 

ii) : El = — i 02 . (2.40) 

Again these two solutions define two lines in the space {E2, 02}- The area between the lines 
corresponds to the parameter cases when the periodic orbit II is hyperbolic due to tangent 
bifurcation. 



D. Symmetry breaking 

The scenaria of the tangent bifurcations of periodic orbits I and II from the last two 
chapters can be studied numerically by use of Poincare maps. For that we always fix the 
total energy of the system, and plot the pair {Xi,Xi} if the conditions X2 = 0, X2 > 
are fulfilled. First we use this circumstance for checking the validity of our results. We 
find excellent agreement. More precisely we test our analytical calculations by evaluating 
numerically the Hill's eigenvalue problem ( p.l9| ) and ( 2.33|) for the two orbits. We find the 



tangent bifurcation, i.e. the merging of the two Floquet multipliers at +1 exactly for the 
predicted parameter combinations. 

An example for the bifurcation of periodic orbit I is shown in Fig. 1 for the parameter 
cases f2 = 1, fs = 3, f4 = 1, 02 = 0.01, 03 = 1, 04 = 1. According to ( |2.31| ) the bifurcation 
energy is given by E^ = 7.407- 10^^. In Fig. la we show the Poincare map of the surrounding 
of the periodic orbit I for Ei = 7 - 10~^ and in Fig. lb for Ei = 7.7-10"^ (note that we did not 
show Poincare maps too close to the bifurcation point only because of the computational 
time which increases to infinity at the bifurcation point because the winding number comes 
very close to unity). Clearly in Fig. la the orbit I is still of elliptic character, whereas in 
Fig. lb the bifurcation already occured, giving rise to the birth of two new periodic orbits 
labelled by A and B in the plot respectively. These bifurcated orbits are separated by a 



separatrix which contains the old periodic orbit I. 

The bifurcation of periodic orbit II is demonstrated in Fig.2 for the parameter cases 
V2 = 1, = 0, Vi = 1, 02 = 0.001, 03 = 1, 04 = 0.1. According to (|239|) , (pliq) the 



bifurcation occurs for E2 = 7.0258 ■ 10 ^. We plot the imaginary part of the Floquet 
multiplier evaluated from equation (|2.33|) as a function of the energy i?2- The obtained 



value for the bifurcation energy is E2 = 7.042 ■ 10^^, which gives a deviation of 0.2% of the 
value from our perturbation approach. 

The considered system of two coupled oscillators exhibits permutational symmetry 
H{Xi,X2) = H{X2,Xi) if we require 02m+i = 0. Still the tangent bifurcation of orbits 
I and II occur. If we define the permutational operator P as 

Pg{X,, X2, Xi, X2) = g{X2, Xi, X2, X,) 

then the normal coordinates are transformed as 

PY, = Y, , PY2 = -Y2 . (2.41) 

Let us consider the bifurcation of orbit I. This orbit corresponds to a closed loop in the 
phase space of the system. Its position is restricted to the subspace {Y2 = F2 = 0}. The 
trajectory of orbit I evolves along the loop parametrically with time. If a tangent bifurcation 
occurs, then the bifurcated orbits will correspond to slightly deformed loops of the original 
one. These deformed loops will move out of the subspace {Y2 = Y2 = 0} as we have 
shown. Let us assume that the deformed loops are invariant under permutation. Applying 
P to such a loop is equivalent to a reflection at the subspace{y2 = 0,1^2 = 0}. It is 
clear that it is impossible to construct such a deformed loop without crossing the subspace 
{Y2 = 0, 1^2 = 0}. But it is forbidden to cross this subspace, since any point from it belongs 
to a periodic orbit I (parametrized by its energy). Thus we conclude that the deformed 
loops can not be invariant under permutation. So the bifurcated orbits break permutational 
symmetry. Then there have to exist at least two bifurcated orbits, which are transformed 
into each other by applying P. Indeed in Fig. lb we observe two bifurcated orbits. We have 
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tested the permutation of these orbits numerically, and found that both orbits break the 
permutational symmetry and are transformed into each other by applying P. Moreover, 
all quasiperiodic motions surrounding any one of the bifurcated orbits and bounded by the 
separatrix also have broken permutational symmetry. The whole corresponding phase space 
regions are connected to each other by P. 

It is straightforward to show that also the tangent bifurcation of orbit II creates bifurcated 
orbits which break the permutational symmetry. The scenaria of both bifurcations are 
identical. This symmetry breaking of the dynamics of two coupled oscillators has been 
observed also numerically by Vakakis and Rand ||20|1 . 

Let us summarize what we have shown so far. We studied the problem of two coupled 
oscillators and derived the dependence of the energy of tangent bifurcation of normal modes 
on the model parameters for small energies. We have shown analytically and demonstrated 
numerically that the new bifurcated orbits will break the permutational symmetry of the 
system (if it existed). 



III. THE CASE OF THE LATTICE 

Let us now turn to the main subject - the bifurcational problem of periodic orbits of a 
lattice. We will consider a one-dimensional lattice with nearest neighbour interaction and 
one degree of freedom per unit cell. We will then show that the derived results apply also to 
lattices with larger interaction ranges. The problem of lattices with higher dimension will 
be discussed. 

The Hamilton function of the lattice is given by 

N 



1=1 



(3.1) 



The potential functions V{z) and ^{z) are again represented in the form (|2.2|) , (|2!^ ). We 
assume periodic boundary conditions 

= Xi^N , Xi = X;+7V ■ 
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Then system exhibits permutational symmetry. The permutational operator P is 

defined by 



Pgi^l: ^2, ^N, ^1, ^2, ■■■,^n) — 9(^2, ^3, -^A^; -^1; ^2, ^3, ■■■y^N, ^l) ■ (3.2) 



Clearly = 1 and PH = H. 

Let us introduce normal coordinates 



N 



1=1 

The wave number q can take any of the values 



2tt 

q = ^n , n = 0,l,2,...,(iV-l) 



The inverse transform of p.3|) is given by 



(3.3) 



(3.4) 



The equations of motion for the normal coordinates Qq read 



N 



1 ini dH 
1=1 



1=1 



AT ^ " dXi 



Using p.lp and ( |2.2| ),( pl3D we obtain the following lengthy expression 



^=3 qi,q2,...,q^-2 



u=l 



TV oo 



fi-i 



E( 



-iq 



l)e-*<?"'Q, 



(3.5) 



(3.6) 



Here ujq abbrevates the eigenfrequencies of the linearized (in Qq) equations of motion and is 
given by the dispersion relation 



= f 2 + 402sin^ Q 



(3.7) 



Let us give the solutions for two periodic orbits of the considered lattice, which correspond 
to the band edge plane waves {q = Q and g7v/2 = tt) in the limit of small energies: 
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I : Qq^o = , Qq=o = - 51 Vl^Qq=0 : (3-8) 

fj.=2 

oo 

II : Q,^^ = , Q,,/, = - Y: ^.^U • (3-9) 

/i=2,4,6,... 

The paramter is given by 

(note the difference to ( |2.9| ). In case II we have to demand (as in the previous chapters for 
the two oscillators) 

II : V2m+1 = 

in order to be able to continue the upper band edge plane wave to finite energies in the 
given form. The terms 02m+i can be in general nonzero, but they simply do not contribute 
to ( p.9|) because of the odd symmetry of the upper band edge plane wave [|15| . 



A. Tangent bifurcation of orbit I 

Let us consider a small perturbation {5q} of the periodic orbit I 

Qq Qq + Sq . 

Linearizing the equations of motion for the perturbation we obtain 

oo 

5q = -COpq - - ^Mq=0Sq • (3.10) 

The analogy to the problem of two coupled oscillators is striking (cf. ( p.7|) ,( pT8D , (|2.19|) ). 
For q = equation ( p. 101) describes the continuation of the periodic orbit itself. All other 
perturbations do not couple with each other, so that we can consider ( p.lOj) for each value 
of q separately. If we increase the energy 

Ej = lQl=o + V{Qq=o) (3.11) 
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(note that the periodic orbit solution has to be inserted in then the first tangent 

bifurcation will occur ii qi = 2t[ /N is choosen in equation ( p.lO| ). For large values of N we 
have 



2 



and consequently obtain (using ( p^.lOD -( P?T3D and ( [^.2(j| )-( pT^ )) the following bifurcation 
energy Ef. 

i) : (^2 = , Ej arbitrary , (3.12) 

• - iv^io^l^g^ • ^^-^^^ 

Since we have to require positive values for Ej and 02 tangent bifurcation can take place 
only if 



f 4 10 V 



< (3.14) 

f 2 9 vi 

in full accordance to equation (|2.32|). Note that this condition is only neccessary, since we 



have also to require f 2 7^ (in the case f 2 = the dependence of the periodic orbit solution 
on Ej becomes different; since this is a nongeneric case, we ignore it). Condition ( p. 14] ) is 



equivalent to the condition, that the frequency of the lower band edge plane wave decreases 
with increasing energy. In other words, the necessary condition for a tangent bifurcation of 
the lower band edge plane wave is the repelling of its frequency from the linear spectrum 



( p.7|) with increasing energy. 

The energy Ej is on the scale of total energy per particle (cf. (|0|),( pl3D ,( p.ll| )). Con- 
sequently the energy threshold of the tangent bifurcation decreases as A^~^ on the scale of 
energy per particle. In the limit N ^ 00 the threshold goes to zero. 

It is important to note, that if all variables Xi are real, the normal coordinates as defined 
are complex. However since the equation ( |3.1CI| ) is linear in 5q and since the mode Qq=o is real, 
we can consider the perturbational problem always separately for the real and imaginary 
parts of a perturbation 5q. 

14 



B. Tangent bifurcation for orbit II 



We again derive the linearized equations for the perturbation 5q 

oo 



OO 

-I 



i E (/^ - l)'^.,...^?^-;/'^^ • (3.15) 

/i=3,5,.-- 



Here we have used the following notations: 

W = ^M + sin^(^)2^<^. , (3.16) 

h,,., = -sin(ga,6)2>^ . (3.17) 

The two wave numbers Qa^b are related to each other by 

gfe = ga ± TT mod27r . (3.18) 

In contrast to the previous cases we have a coupling between a pair of normal coordinates 

( 3. IS ) in (p. 15 ). Notice that the coupling term is given by the second sum on the right hand 



side of (|3.15|) and is zero if 02m+i = 0. Also this coupling term is proportional to i, which 



causes a mixing of real and imaginary parts of the perturbations. 

Interestingly for the pair = vr, g;, = we obtain again no coupling, because (|3.17|) 



vanishes for both wavenumbers. Consequently for g = tt ( |3.15|) describes the continuation 
of the periodic orbit II. 

1. The case (p2m+i = 

If we assume 02m+i = then the equations for the perturbations 6g^ ^ decouple. In 
analogy to the case of orbit I we have the first tangent bifurcation of orbit II if we consider 
the perturbation g = 7r(l — 2/N). Consequently we obtain for the bifurcation energy Ejj 

i) : 02 = , E'ij arbitrary , (3.19) 
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Again we observe that a bifurcation will take place only if is positive, i.e. only if the 
frequency of the periodic orbit II is repelled from the linear spectrum with increasing energy. 

We can compare these findings with the calculation of Budinsky and Bountis |I6[ , who 
have calculated the bifurcation energy for v^=o, 02 = 04 = 1, = otherwise. They obtain 

= 3.226/iV2 (equation (2.22) from Pf). Our resuh gives = 7i^/{3N^) ^ 3.29/N^. 



Budinsky and Bountis have roughly estimated the eigenvalue spectrum of the corresponding 
Hill's matrix in the original coordinates. We think that this circumstance is the reason for 
the small deviation of their approximate result from our exact one. 

Note that Budinsky and Bountis are looking for the first bifurcation to occur, either 
tangent or period doubling. It can happen that a period doupling bifurcation occurs prior 
the tangent one. In that case only the result for the period doubling bifurcation would be 
given in their paper. In the quoted comparison indeed the tangent bifurcation takes place 
first, thus the apparent coincidence in the results. However for other cases period doubling 
bifurcations can occur first. Then our results can not be compared to the ones of Budinsky 
and Bountis (cf. also the discussion section). 

2. The case 4>2m+i / 

Now the coupling between two perturbations has to be taken into account. We use the 
fact that Qqj^/^ is real and the equations for the complex perturbations are linear. By taking 
the complex conjugate of any of the two equations defined by ( |3.15|) and considering say 



only the real part of 5q^ and the imaginary part of 6q^^ we find in the limit of small energies 
En 

'L,, = + ^^4,qa,tQl^^,)K, - h,t,aQlM/2Ka ' (3-21) 

Note that we have changed the relation between and qb to 

Qb = \qa - 7r| , Qa mod7r/2 . 
The lowest tangent bifurcation energy of orbit II appears for 
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Then it follows 



k,q^ = = -1603 



N 

One could expect that in the limit of large the coupling between and qh vanishes. That 
is indeed so if we require f 2 7^ (linear spectrum is optical-like) but turns out to be wrong 
for the case f 2 = (linear spectrum is acoustic-like). The details of these subtleties are 
given in Appendix A. Here we proceed to the final result for the tangent bifurcation energy 

Eh-- 

i) : Eh = ^ (^2 + 402) (02 + ^^^2) 02 , (3.22) 

-WWI^)^- -2^V (3.23) 
i67i£ 4:1 „ - n 

302K+1604)-64<A^ ^ 

As we can see solution i) (|3.22|) is always positive (since 02 > and V2 > 0) but is not 
dependent on A^. Since we applied perturbation theory, (|3.22|) is correct only in the limit 



of small energies. For larger energies corrections apply. Solution ii) ( p.23| ) is the one which 
gives arbitrarily small bifurcation energies for sufficiently large A^. Since 03 does not enter 
the energy dependence of periodic orbit 11, we then again obtain as the neccessary condition 
for the existence of the bifurcation, that the frequency of the upper band edge plane wave 
has to be repelled from the linear spectrum with increasing energy. However for the case 
f 2 = a more restrictive condition is obtained by demanding 

302(t;4 + 1604) > 6402 _ (^3 24) 

Sandusky and Page |]T3| have considered the case 02 7^ 0, 03 7^ 0,04 7^ 0,0^ = otherwise. 
However in contrast to Budinsky and Bountis no clear finite size studies are done in [ll5[ . 



Instead these authors use dimensionless parameters, which however contain the amplitude 



(energy). Still we can consider equation (24) from [|T5| and by performing the limit of small 
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amplitudes we recover our neccessary condition ( p.24|) . The figures 8 and 9 from |]T5| are 



sliowing a crossover from tlie parameter region wliere ( |3.24| ) liolds into the region where it is 
not vahd. From our analysis it follows that this corresponds to a change from ( p.23| ) to ( |3.22|) 



and consequently to a change from infinitesimally small energy (amplitude) thresholds to 
finite thresholds. Indeed finite amplitude thresholds for the instability have been observed 
in |]l5l after leaving the parameter region where ( p.24|) is valid. 



It is interesting to note that condition (|3.24|) has been obtained with the help of multiple 



scale expansions already in 1972 by Tsurui (equation (4.8), uJl = 4) and more recently 
by Flytzanis, Pnevmatikos and Remoissenet p2| (equation (5.15)) for systems with = 0. 



C. Symmetry breaking 



Using (0),(|0) it follows 



P"Q, = e^'?"g, , q = ^rn , m = 0, 1, 2, (iV - 1) . (3.25) 

First we notice that the permutational operator does not mix the space of the normal 
coordinates. In other words, every subspace {Qq, Qq} is invariant under ( p. 21) . However since 
the phase space is higher dimensional compared to the case of two coupled oscillators, the new 
bifurcated orbits can not be confined to some subspace (in order to use the argumentation 
of the two coupled oscillators case). 

Let us consider the bifurcation of orbit I. We will project a new bifurcated orbit loop 
into the subspace {Qo,Qo,Qi,Qi}- We can use a reference picture in a three-dimensional 
space, where the Z-axis contains both variables Qq and Qq. Now the projected loop can 
cross the Z-axis (in contrast to the case of two coupled oscillators), because at the moment 
of the crossing the original loop can still not intersect with the subspace {Qo, Qo}- However 
we notice that every loop (also the projected one) have a certain direction of rotation of the 
trajectory. If the loop is invariant under P, so is its projection. The direction of rotation 
can either be conserved or change sign under permutation. So if we choose n = N/2, n even, 
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in (|3.25|) then the direction of rotation is conserved. Then we observe, that the constructed 



operation P^/^ leaves the periodic orbit I, but switches the sign of the coordinates of the 
relevant perturbation Qi. Applying the permutation to the projected loop is equivalent to 
reflecting the loop at the subspace {Qo,Qo} (the Z-axis in the reference picture). Clearly 
this reflection switches the direction of rotation. Consequently our initial assumption was 
wrong, and the bifurcated orbits indeed break permutational symmetry. It is straightforward 
to show that also the new bifurcated orbits in the case of a tangent bifurcation of orbit II 
break permutational symmetry. 

Then in general the new bifurcated orbits will be not invariant under P, and thus the 
bifurcated orbits will correspond to band edge modes, which are spatially modulated by the 
relevant perturbation. This perturbation can be represented as a wave with a wavelength 
equal to the length of the system (see Fig.3). Note that the modulation referrs to the additive 
perturbation, so the bifurcated orbits do not vanish in real space for certain coordinates Xi. 
They are thus formally similar to standing waves, but also inherently different. Applying 
a permutational operation, we will simply shift the modulated object through the lattice. 
Consequently we arrive at the result, that after the tangent bifurcation of a band edge plane 
wave (orbit I or II) takes place, at least new periodic orbits bifurcate from the original 
one. 

It is much harder to argue about the exponential character of spatial localization of the 
bifurcated orbits. For any large but finite A^ the qualitative form of the bifurcated orbit 
as shown in Fig.3 is correct if the amplitude of the perturbation is small compared to the 
amplitude of the bifurcating orbit (I or II) at the bifurcation. The frequency of the orbits I 
(II) and of the newly bifurcated ones are outside the linear spectrum. Let us consider the 
analytical continuation of a newly bifurcated orbit. Its frequency will follow (with deviations 
of higher order) the one of the original periodic orbit (I or II), i.e. A(u;^) ~ E/N. So we can 
assume that the frequency is further repelled from the linear spectrum. The amplitude of 
the original periodic orbit grows like (E/NY^"^. The amplitude of the relevant perturbation 
however grows much faster. That is due to the fact, that at the bifurcation at least A^ new 
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periodic orbits occur. So in terms of catastrophe theory we would be confronted with an An 
singularity |^3[. The amplitudes of the relevant perturbations will grow as {E/N — E'^)°', 
where a < 1/2. Thus the amplitude of the relevant perturbation (describing the deviation 
of the new periodic orbits from the original one) will grow much faster than the amplitude 
of the original orbit as a function of energy. Consequently we can expect, that the overall 
amplitudes in region II in Fig. 3. will be drastically lowered, whereas in regions I and III 
they will grow with energy. This gives us a tendency towards localization, but not more 
than that. If it indeed localizes more with increasing energy then the spatial decay will be 
exponential in region II, because the frequency of the bifurcated orbit is outside the linear 
spectrum. Clearly it is necessary to explicitely take into account all nonlinearities in order 
to give a more precise answer on the question of localization. 



IV. DISCUSSION OF THE RESULTS 

We have studied the tangent bifurcation properties of band edge plane waves in the limit 
of large system size. We have restricted the consideration to one-dimensional lattices with 
nearest neighbour interaction and one degree of freedom per unit cell. Clearly the results 
can be generalized. 

First we consider the case of higher dimensional lattices. Breathers have been numerically 



observed in [§, pj, p5| and analytically obtained in |jTj| for selected systems. Again the 
used approach can be applied. It is essential to know the properties of the linear spectrum. 
Then one has to evaluate the system size dependence of the relevant perturbations (the ones 
with the eigenvalues closest to the band edge frequency) and to repeat the calculations. 
Without performing this calculation (which is anyway easy to do using the results of this 
work) we mention that the system size N of the one-dimensional lattice (i.e. the number of 
degrees of freedom) has to be changed to N^^'^ in the expressions for the bifurcation energy 
derived in this paper, where d is the dimensionality of the lattice. So except for some changes 
in the formula no qualitative differences appear. Still there is a problem - as one can see, 
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the total bifurcation energy is independent of for d = 2 and even increasing with N^^^ 
for d = 3. Using [T^ we know that in weakly coupled anharmonic oscillator arrays discrete 
breathers exist at arbitrary small energy for zero coupling, and thus at small energy for weak 
coupling. Consequently the bifurcation of the band edge plane wave for d = 2,3 might be 
not directly related to these solutions. 

We can only suspect, that other channels of breather birth through bifurcations exist. In 
any of the considered systems standing waves occur in the linearized case which correspond 
to resonant tori (two identical frequencies). These standing waves already in the linear case 
break the permutational symmetry. The interesting cases correspond to wavenumbers closest 
to the band edge plane wave. The standing wave would then have two distinct knots along 
either spatial direction, since the original waves had a wavelength equal to the system size. 
These resonant tori will most probably transform into chains of regular islands separated 
by a separatrix, once the nonlinearities are included. The newly bifurcated orbits can shift 
their frequencies out of the linear spectrum with increase of energy. Since they were close to 
the band edge plane wave it is likely, that the energy dependence of the frequencies of these 
new orbits is similar to the one of the band edge plane wave. If this scenario is true, then 
one could still predict the existence of breathers in higher dimensional lattices by analyzing 
the energy dependence of the band edge plane wave frequency. A logical question would 
be, whether the same effect of breather birth through resonant tori (rather then through 
tangent bifurcation of a band edge plane wave) takes place in the d = 1 systems, which were 
considered in the present paper. At the present the answer is, that it might well be. Still 
there exist breathers at infinitely small energies which are bifurcating from the band edge 



plane wave, as shown in this work, and as also observed in numerical studies [0 . One could 
argue that for the A^ = 40 example in Fig. 8 of the dissapearance of the discrete breather 
at the frequency of the band edge plane wave (of same energy) could be nearly reproduced 
by the breather merging with a standing wave which has a difference of about 1/1600 in 
frequency compared to the band edge one. But the main problem with the standing wave 
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is, that the bifurcating objects will have the symmetry of two breathers. At the present we 
do not see another channel for a single breather birth via bifurcation rather than the band 
edge plane wave bifurcation discussed in the present paper. 

Secondly we consider the problem of larger interaction range. Again changes in the 
linear spectrum have to be considered, and for the upper band edge plane wave (orbit II) 
also some changes in the derivation have to be expected. For the case of weak next-to 
nearest neighbour interaction (compared to the nearest neighbour one) a smooth change 
of the derived results will apply. Especially we observe, that the tendency is to increase 
the distances between the eigenvalues of the linearized spectrum (at fixed A^) and thus 
to increase the bifurcation energy. Thus increasing the interaction range corresponds to a 
gradual suppression of bifurcation and localization. That seems to be plausible, since in 
certain limits of mean-field type interaction no localization should occur. 

The analysis of systems with several degrees of freedom per unit cell seems to be more 
complicated, allthough the presented approach can be used. It becomes a matter of right 
dealing with polarization vectors etc. Here we can only expect that tangent bifurcations of 
band edge modes will appear if the neccessary condition of the mode frequency repelling 
from the linear spectrum (band) occur. 

In this work we have studied only tangent bifurcations of band edge plane waves. These 
cases have to be clearly separated from e.g. period doubling bifurcations (cf. ||15|). Espe- 
cially the subtle connection between tangent bifurcation of band edge plane waves and the 
appearance of new periodic orbits which break permutational symmetry and exponentially 
localize in the configurational space is not easily transferred to the case of period doubling 
bifurcations. Since period doubling bifurcations will occur only for the upper band edge 
plane wave, the bifurcated orbits will have frequencies located inside the linear spectrum. 
Thus we doubt that any tendency for exponential localization can occur in these cases. 

Let us estimate the bifurcation energy of a period doubling bifurcation of periodic orbit II. 
The relevant perturbation is given by g = 7r/6 for the one- dimensional case (the upper band 
edge of the linear spectrum is given by 2\f^2-, the relevant frequency of the perturbation has 
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to be close to \f^2) ■ If ^ 1 6 is an integer, then we have to expect a period doubhng bifurcation 
at zero energies, since the normal modes resonate already in the linear problem. If is 
noninteger, then the distance of the closest eigenvalue of the linear spectrum to the value 02 
is of the order Note that this is different from the tangent bifurcation cases considered 
above, where the distance of the closest eigenvalue to the squared band edge was always of 
the order 1/N'^. Thus we conclude that for the generic cases period doubling bifurcations 
occur at energies E'' ~ (cf. Budinsky and Bountis [Q, equation (3.8)) whereas the 
tangent bifurcation occurs at energies E'^ ~ l/N"^. So generically tangent bifurcations 
appear at lower energies than period doubling bifurcations. Then we can conclude, that 
with increase of energy first exponentially localized orbits (discrete breathers) are generated 
through tangent bifurcations. 

If we consider a nonlinear Klein-Gordon lattice (optical-like linear spectrum, V2 7^ 0) 
then for small bandwidth (02 '^^2) no period-doubling bifurcations occur at small energies. 
Since generically any nonlinearity will repell the frequency of one of the two band edge plane 
waves with increase of energy, tangent bifurcation always occurs. This result is in agreement 
with MacKay and Aubry |T^, who have shown that discrete breathers exist for these lattices 
for any nonlinearity, provided the bandwidth is small. 

Let us give an example where no discrete breathers are expected. The one- dimensional 
To da chain is characterized by Vfj_ = 0, 02 = 1? 03 = ^1/2, 04 = 1/6. Clearly (|3.24|) 



is not satisfied (note that in that case the frequency of the upper band edge plane wave 
is repelled from the linear spectrum with increase of energy, so the repelling condition is 
a necessary condition, but not a sufficient one). Thus in the limit of large system size no 
tangent bifurcation occurs at small energies, and no discrete breathers should occur. This 
circumstance should make clear, that discussions about 'nonexistence of breathers in systems 
with realistic potentials' are simply based on the fact, that a potential called realistic can 
be choosen in a way that (|3.24|) does not allow for tangent bifurcation of an upper acoustic 
band edge phonon. 

A few words should be added in order to explain the subtleties of the result of the previous 
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paragraph. If we consider a system with an acoustic-type hnear spectrum (hke that of the 
Toda chain) then we deal with equations ( |3.22| ) and ( |3.23| ). The first tangent bifurcation 



occurs at the value of energy which corresponds to the smallest nonnegative outcome of any 
of both equations. If ( |3.24| ) is fulfilled, then ( |3.23| ) gives positive energies, which become 



arbitrarily small as the size of the system is increased to infinity. Consequently the new 
bifurcating orbits have zero amplitudes and energy densities outside their modulation center 
and eventually evolve into discrete breathers with increase in energy. If however ( p. 241) does 



not hold (e.g. as in the case of the Toda chain) then ( p. 231) gives negative bifurcation energies. 
So we have to use ( p.22| ), which gives finite energy values even for infinite system size. Despite 
the fact that ( |3.22|) is only a result of perturbation theory, and corrections apply for finite 



but not too small energies, we can conclude for sure that the tangent bifurcation energy in 
this case stays finite with increase of the size of the system. So the new bifurcated orbits will 
have finite nonzero amplitudes and energy densities outside their modulation center even if 
being infinitely close to the bifurcation. Those objects can not be called discrete breathers, 
because their amplitudes can not decay exponentially in space, and because their total energy 
on an infinite lattice will be infinite. Rather they correspond to discrete breathers 'sitting' 
on a nondecaying carrier wave. We suspect that precisely those solutions were observed 
in Figs. 8, 9 in [|l^ but were not identified as being different from discrete breathers with 
exponential spatial decay. So the result of the previous paragraph is that there exist no 
discrete breathers in a Toda chain. But the next result and prediction is that there exist 
discrete breathers on carrier waves in the Toda chain (this last statement should be still 
considered to have some uncertainty, because of the mentioned use of perturbation theory 
for deriving ( p.22|) ). 

Nonlinearity shifts frequencies of periodic orbits with change of energy. For an optical 
band no matter whether the shift results in an increase or decrease of the frequency, tangent 
bifurcation will take place nearly always - either for the upper or for the lower band edge 
phonon. In the case of an acoustic band we need increase of frequencies in order to obtain 
tangent bifurcation, and even this condition is only a necessary one (cf. the two previous 
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paragraphs). Thus acoustic band related breathers will occur rarely compared to the optical 
band related ones. This circumstance is in accordance with numerical findings p8[] 
where periodic modulation of masses (diatomic chain instead of monoatomic chain) allows 
for tangent bifurcation of band edge plane waves of the optical band (which is absent in the 
monoatomic case, as were the breathers for the choosen parameters). 

From our results we can follow, that breathers will be likely to occur in lattices with 
optical bands in the linear spectrum, which are narrow and well separated from the rest 
of the spectrum (we have recently learned about a numerical study of Bonart, Meyer and 
Schroder p9|], where similar guides have been successfully applied in order to obtain discrete 
breathers on surfaces of 3d crystals). Of course the nonlinearity has to be strong enough - 
that can be accomplished either by not too low temperatures (especially for systems with 
structural phase transitions) or by local excitation of the lattice using e.g. laser pulses, 
so that sufficiently large amplitudes occur, which should then provide for strong enough 
nonlinearities in the forces. Breathers will certainly infiuence the thermal conductivity 
acting as strong scatterers of small- amplitude phonons [^] and may also infiuence other 
transport properties of other degrees of freedom (e.g. electrons) if the coupling is strong 
enough. 

Finally we want to note, that for systems with acoustic-like linear spectra (say = 
0) and 02m+i 7^ numerical simulations of discrete breathers in one-dimensional chains 
always lead to the presence of static displacements (dc terms) in the solution Xi{t) which 
seemingly decay linearly in space, due to the periodic boundary conditions |15[, |^ 



In the limit of large system sizes the slope goes as and thus vanishes in the limit 
N ^ oo. This result can be immediately obtained using the fact that the decay rate of the 
dc term is zero for an infinite lattice with an acoustic linear spectrum [|10|. Thus discrete 
breathers in the mentioned lattice types are similar to a kink+localized vibration for the one- 
dimensional case, and to a kink bubble+localized vibration in higher dimensional lattices. 
This statement steems from the observation, that far away from the center of the breather 
the static displacements can be described as a strain field caused by a point-like source 
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(in full analogy to the corresponding Maxwell equation for the electric field) ||3^. Then 
the strain field will be radial and isotropic, and its decay can be calculated using Gauss's 
theorem. The resulting decay is ~ l/r'^"^, where d is the dimension of the system. So for 
three-dimensional lattices we obtain decay, for two-dimensional systems 1/r decay and 
for one- dimensional systems non-decaying strains. Consequently even in three-dimensional 
lattices these discrete breathers (which bifurcated from an acoustic band in the presence of 
02m+i 7^ 0) are composed out of a dynamical (time-periodic) solution decaying exponentially 
fast, and a strain field decaying as If we choose ^2 7^ (i.e. if we consider a breather 

which bifurcated from an optical band) then the spatial decay of the dc terms becomes 
exponential. Similar to that we find that the tangent bifurcation energy of the upper band 
edge plane wave behaves also differently for the two cases ( |3.23| ). Performing first the limit 



— > 00 we thus find that the two different classes of lattices (acoustic-like or optic-like 
linear spectra) have no simple contact. This is connected to the circumstance, that the case 
of acoustic-like linear spectra implies the existence of an additional integral of motion - the 
total mechanical momentum. In the case of systems with more than one degrees of freedom 
per unit cell it can happen, that tangent bifurcation appears for an optical band edge plane 
wave. However the nonlinearities can couple the band variables of the optical band with the 
acoustic band variables. Then the existence of a local strain due to the discrete breather 
will still cause a nonexponential decay of the strain as described above. Thus we agree with 
the conclusion in (where a nondecaying strain was found in a one- dimensional diatomic 
lattice - with the breather bifurcating from the optical band) that the existence of breathers 
in crystals should be in general accompanied with anomalies in the thermal expansion. 

There exist approximate methods of testing whether a lattice allows for the existence 
of breathers 0, 0, [^]. These methods use the hypothetical existence of breathers and 
construct then energy dependencies of their frequencies. If the breather frequencies are 
attracted by the linear spectrum, then the conclusion was that breathers most probably do 
not exist. These predictions were usually successful. However this method of argumentation 
is not well-defined, since it uses objects which are not defined (as long as we do not know 
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whether breathers exist, we can not calculate any property of a breather). In the present 
work instead a well-defined and clear criterion is given - whether or not the frequency 
of a band edge plane wave is repelled from the linear spectrum with increase in energy. 
There should be even experimental methods of testing this property of a lattice, if the first 
corrections to the linear spectrum due to nonlinearity can not be well estimated theoretically. 



V. SUMMARY 

We have shown, that band edge plane waves of nonlinear Hamiltonian lattices undergo 
tangent bifurcations. A necessary condition for that is the repelling of the frequency of the 
plane wave from the linear spectrum with increase of energy. If this condition is fulfilled, 
then the bifurcation energy (in units of one-particle energy) scales as N^'^/'^ where d is the 
dimensionality of the lattice. The new bifurcated periodic orbits break the permutational 
symmetry of the lattice (if periodic boundary conditions are used). The shape of the new 
bifurcated orbits corresponds to a localized vibration. It is argued that the spatial decay 
of the amplitudes becomes exponential with further increase of energy. Thus we confirm 
the numerical findings of [jl5l , who have shown that the new bifurcated orbits are related to 



discrete breathers in finite chains. In accord with the strong evidence that discrete breathers 
are generic solutions of nonlinear lattices we find that tangent bifurcation of band edge plane 
waves is generic too. 
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APPENDIX A: PERTURBATION THEORY FOR BIFURCATION OF ORBIT II 

IN THE LATTICE 



Here we want to scetch the derivation of the critical hnes of tangent bifurcation for 
periodic orbit II (upper band edge plane wave) for 02m+i 7^ 0. We start from the set of two 
equations (|3.21|) and rewrite them in the following way: 



X 



nl + 2F^cos{2ujt) X - 2Aacos{ujt)y , 
nl + 2F^cos{2ujt)] y - 2Aacos{Lut)x . 



(Al) 
(A2) 



Here we used the different notations 



X = ^qa , y = ^> 



a = (j): 



The other parameters are defined in the limit of small energies of the band edge plane wave 
as 



2 2 



A = aiE 



1/2 



V2 + 



2V2 

47r2 , 



2 



00. 



y 



f 2 + 



47r2 



(A3) 
(A4) 

(A5) 
(A6) 

(A7) 
(A8) 



We expand the solutions of the differential equations (|Al|) , ([A^) into series of powers of E^^'^: 



x{t) = xo{t) + E^l'^xi{t) + Ex2{t) + ... , 
y{t) = yo{t) + E'l''y^{t) + Ey^it) + ... , 
nl=uj'' + giE^'^ + g2E + ... . 



(A9) 
(AlO) 
(All) 



Note that we did not expand Qy, since in the given problem this frequency is far away from 
the frequency of the periodic orbit u. 
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Inserting (|A9|) - (|A11|) into (^),(^) and sorting with respect to powers of E^/^ we obtain 
in lowest order 

io = -t^^xo , y'o = -ujyVo . (A12) 
Since we are looking for solutions periodic with 2^/00 it follows 

xq = acos{ujt) + bsm{u!t) , yo = . (A13) 

In order E^^'^ we get 

x\ = uP'xx — giXo , y'l = —^yVi — 2aaicos{ujt)xo ■ (A14) 
The solution reads 

^1 = 0, xi = , (A15) 

yi = K + ccos{2ujt) + dsm{2ut) , (A16) 

aai c d aai / a i ^\ 

^ = ' ~ = I = 2 ■ 

ujy a Auj^ — ujy 

Note that the solution xi = is not strictly required at this level. However we would have to 
do so at the next level, so we remove it already here. The constant n is inverse proportional 
to Uy. This circumstance will lead to the difference in the final results for the optical like 
spectrum {Uy finite) and the acoustic like spectrum {Uy ~ 
In order E we obtain 

X2 = —uj'^X2 — g2Xo — /2,xCos(2co't)xo — 2aaicos{ujt)yi , (A18) 

i/2 = -^2Z/2 . (A19) 

Equation ( [A19| ) simply requires y2 = 0. Removing the secular terms from equation ( [A18| ) it 
follows 

(A20) 

ii) ■ E'^ = %<P^ ^ X • (A21) 



y 
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Note that o? ~ l/N"^. Thus we obtain that the A^-dependence is removed from equation 
( |A2CI| ). Now we see the cause for the mentioned difference between the optical and acoustic 
spectra. If uj"^ stays finite for large A^, then the whole second term in the denominator of 
the right hand side of equation ( [A21| ) scales to zero. In the case of the acoustic spectrum 
~ and the same term in ( |A21D gives a finite contribution in the limit of large N . 
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FIGURE CAPTIONS 



Fig.l 



Poincare plot for the system of two coupled oscillators. The condition is X2 = 0, > 0. 
The parameters are given in the text. Only a part of the available phase space is shown, 
which includes the relevant periodic orbit I. 

a) The energy is below the bifurcation; 

b) The energy is above the bifurcation. A and B labels the two new periodic orbits, which 
appeared due to the bifurcation of the periodic orbit I. 



Fig.2 



The imaginary part of the Floquet multipher (absolute value) of the Hill's equation as a 
function of the energy E2 for periodic orbit II and parameters as given in the text. At 
energies larger than the bifurcation energy the imaginary part vanishes. 



Fig.3 

Schematic representation of the relevant perturbation of a periodic orbit I due to the sym- 
metry breaking. The x-axis represents the normalized spatial variable l/N , and the y-axis 
the amplitude of the solution in arbitrary units. The dashed line indicates the amplitude 
distribution of the original periodic orbit I. 
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